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Abstract 

It is shown that in polar geometry and normal incidence the 2x2 matrix 
technique - as discussed in detail in a preceeding paper [1] - accounts correctly 
for multiple reflections and optical interferences, and reduces only in the case of 
a periodic sequence of identical layers to the Fresnel formula of reflectivity which 
in turn is the theoretical basis of the two-media approach, widely used in the 
literature to compute magneto-optical Kerr spectra. As a numerical example 
ab-initio calculations of the optical constants for an fee Pt semi-infinite bulk 
using the spin-polarized relativistic screened Korringa-Kohn-Rostoker method 



show very good agreement with experimental data. 
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INTRODUCTION 



In the last few years magneto-optics started to become of prime interest in dealing with mag- 
netic multilayer systems and magnetic nanostructures (submonolayer coverages of substrates 
with magnetic material): Kerr measurements not only turned out to be one of the standard 
experimental tools applied, new applications, in particular in the context of time-resolved tech- 
niques, are presently believed to lead to fast magnetization switching devices. In many cases, 
however, theoretical descriptions even for time-integrated magneto-optical effects are lagging 
behind the experimental efforts mainly because relevant schemes to deal with semi-infinite sys- 
tems - not to speak of nanostructured materials - are not commonly used. Bulk-like approaches 
(assumed three-dimensional periodicity) and the so-called two-media model [4] (assumed ho- 
mogeneity) are still considered to be sufficient to deal with magnetically inhomogeneous layered 
systems. 

In a recent paper, [1] the authors showed that only by including multiple reflections and 
optical interferences, e.g., via the 2x2 matrix technique [2,3], realistic ab-initio magneto- 
optical Kerr spectra for semi-infinite layered systems can be obtained. It is the purpose of 
the present paper to prove first analytically and then numerically that only in the case of 
periodic layered systems, i.e., in by definition homogeneous systems, the 2x2 matrix technique 
reduces to the two-media approach. Formulated oppositely, this implies that the two-media 
approach is strictly valid only for this kind of systems. This is illustrated for fee Pt viewed as a 
periodic layered system, because Pt frequently serves as substrate for magneto-optically active 
multilayers. [6,7] 

THEORETICAL APPROACH 

— * 

Assuming polar geometry and the magnetization M p in all layers p to point along the z- 
direction, the layer-resolved permittivity in cubic, hexagonal or tetragonal systems is given 
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provided that a possible anisotropy in the diagonal elements can be neglected, i.e., e p z ~ e p x . 
In case of normal incidence, n px = h py = 0, four electromagnetic beams corresponding to the 
complex refractive indices 



~(3) 



n 



-n (1) 

^(2) 



pz 



-n 



pz 



propagate in layer p. By considering harmonic fields 



A(z, t) = „4exp [i (qz — cut)] exp (St) = .4* exp [i (qz — cot)] 



(2) 



such that 5 > describes the interaction between the layered system and its neighborhood, 
q = q$n is the complex wave vector with q being the propagation constant in vacuum, and 
u = to — id, beams 1 and 2 propagate along — z, beams 3 and 4 along +z. The surface reflectivity 
matrix, see in particular the appendix of Ref. [1], 
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which relates the incident electric field to the reflected one, 

n- 1 
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-Rsurf = (B vac + V N A *) (-Bvac ~ T> N A *) , 
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is obtained as defined in Eq. (37) of Ref. [1], in terms of reflectivity matrices TZ p: see Eq. (32) 
of Ref. [1], 



K p = (B p + Pp_0 1 (B p - Pp_0 , p = 1, . . . , N , 
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and propagation matrices C p , see Eq. (30) of Ref. [1], 
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where d p refers to the thickness of layer p. According to Eqs. (28), (32) and (35) of Ref. [1] the 
2x2 matrices occurring in Eq. (4) are defined as follows 



v p = b p [c p - (Cp)' 1 n p ] [c p + (c p y l n p ] 1 , 
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1Z 7v results then recursively starting from the substrate reflectivity matrix 1Z = 0. It has 
been shown in quite some detail in the appendix of Ref. [1] that the layer-resolved reflectivity 
matrices 1Z P are all diagonal, 
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where r p ± is the complex reflectivity coefficient of the right- and left-handed light in layer p. 
This implies that the right- and left-handed circularly polarized components of the incident 
linearly polarized light, once they arrived at the surface layer, propagate independently within 
the system such that after the first reflection, they become immediately elliptically polarized. 



SPECIAL CASE OF HOMOGENEOUS SYSTEMS 

In principle in a finite periodic layered system (N can be large but finite; for matters 
of simplicity a simple parent lattice [5] is assumed) all layers have identical layer-resolved 
permittivities. This implies in turn that also all matrices B p in Eq. (8) are identical, i.e., are of 
the form 
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with n± = \fs±. The recursion relation in Eq. (5) reduces therefore to 



ll p = (l+ Q p -i) 1 (J - Qp-x) , p = 1, 
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where X denotes the 2x2 identity matrix and 

Qp-i = [Cp-i - {Cp-i)' 1 TZp-i] [Cp-i + (Cp-i) 1 TZ p -i\ . (11) 

This finite periodic layered structure has to be properly matched to a semi-infinite system 
(substrate) of the same material. Inserting TZ = into Eq. (11) for p = 1, yields Q = T, which 
substituted into Eq. (10), immediately proves that also TZi = 0, and so on. In conclusion, for 
a periodic layered system, 

7Z P = 0, Vp=l,...,N. (12) 

^From an optical point of view a periodic layered system behaves like a system with no 
boundaries. Viewed oppositely, Eq. (12) shows that there is a boundary in-between two adjacent 
layers, if and only if, the respective layer-resolved permittivities differ. 

^From IZn = follows immediately that V N = B with B being defined in Eq. (9); hence 
Eq. (4) directly yields the surface reflectivity matrix 
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i.e., r xx and r xy in Eq. (3) assume the following values 

n + n_ — 1 _ . h + — fL 



(n+ + l)(n_ + l) ' rxy \h + + 1) (n_ + 1) ' 

In the case of periodic layered systems the complex reflectivity coefficient of the right- and 
left-handed polarized light is therefore given by 

n±-l M A , 

r± = r xx =F ir xy = — , (14) 

n± + l 

a relation, which is known in the literature as the Fresnel formula for s -polarization and normal 
incidence. [9] Eq. (14) then leads directly to the well-known formula for the complex Kerr angle 
in the two-media approach, [4] 

where 6k is the Kerr rotation angle and 6k the Kerr ellipticity. 



Clearly enough in the case of (homogeneous) periodic layered systems the 2x2 matrix 
technique and the two-media approach provide identical Kerr spectra. However, by using 
the two-media approach for calculating Kerr spectra of inhomogeneous layered systems, such 
spectra are governed almost exclusively by contributions from the surface layer. 



APPLICATION TO FCC PT 

In Figs. 1-3 the optical constants of fee Pt bulk as calculated via the 2x2 matrix 
technique for normal incidence and different surface orientations are compared with available 
experimental data (Ref. [10] and references cited therein). Because fee Pt is paramagnetic, 
the right- and left-handed complex reflectivity coefficients are equal, r + (uj) = f_ (uj) = f(uj), 
and the reflectance R(uj) given in Figs. 1 is in fact f (uj) 2 . In accordance with Eq. (2) the 
refractive index n (uj) and the absorption coefficient k (uj) in Fig. 2 simply correspond to the 
real and imaginary part of the complex refractive index n (uj) = yje (uj) = n (uj) + i k (uj). The 
total (complex) permittivity e (uj) displayed in Fig. 3 is evaluated assuming a periodic layered 
system consisting of N Pt layers on top of a semi-infinite Pt substrate, i.e., 

£ »4e ??(w) ' (i5) 
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with a vq (uj) being the inter- and intra-layer contributions to the optical conductivity and 
uj = uj — id. Eqs. (15) and (16) follow directly from the Fourier-transformed (linear) material 
equations [11] by assuming that each layer p can be viewed as a homogeneous medium and 
therefore the layer-resolved optical conductivity a p (uj) is related to the layer-resolved permit- 
tivity by 

4-7TV 

~ £ p( u ) = l + — d p {uj) ) p , g= i,...,iv. (17) 

UJ 

The a pq (uj) in Eq. (16) are calculated in terms of Luttinger's formula [12] using the spin- 
polarized relativistic screened Korringa-Kohn-Rostoker method (SKKR) for layered systems 
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[13,16-18], contour integration techniques, [8,14] the Konrod quadrature and the cumulative 
special points method for the occurring Brillouin zone integrals. [19] It should be noted that 
Luttinger's formula uses a vector potential description of the electric field and has several 
advantages over the well-known Kubo formula [15]: both, the absorptive and dissipative parts 
of the conductivity tensor are included without using Kramers-Kronig relations, and so are all 
inter- and intra-band contributions avoiding thus a phenomenological Drude term in order to 
mimic the latter contributions. 

In the present study of fee Pt bulk (lattice parameter 7.4137 a.u.) all complex energy and 
/c-space integrals are performed with an accuracy of 10~ 3 (in atomic units) for all surface 
orientations considered. The electronic temperature amounts to T — 300 K. Furthermore, a 
lifetime broadening of 0.048 Ryd ~ 0.653 eV is used and 2 and 37 Matsubara poles in the lower 
and upper semi-plane, respectively, are taken into account. 

As can be seen from Fig. 1 the reflectance R (ui) of Pt as calculated for normal incidence 
and different surface orientations agrees very well with the experimental data. Unfortunately 
the experimental data [10] are not specified with respect to surface orientations, i.e., no surface 
normal is listed. The dependence of all the calculated optical constants on the surface orien- 
tation reflects, however, merely the fact that in principle in spectroscopical experiments only 
physical properties of semi-infinite systems, i.e., of solid systems with surface, are recorded, see 
also Figs. 2 and 3. 

A close inspection of Fig. 1 in the vicinity of very low photon energies reveals that the 
fine structure around 0.5 eV seen in the experimental reflectance is missing in the theoretical 
curves although the maximum in the refractive index as a function of the photon energy is well 
reproduced. This discrepancy between theory and experiments is caused by an intrinsic feature 
of the applied approach: the finite lifetime broadening 5 that enters Luttinger's formula smears 
out possible fine structures for uj < 5 (the value of uj = S is indicated in Figs. 2 and 3 as a 
dotted line). 

In Figs. 1-3 the theoretical data refer to a total of A = 12 Pt layers, because from 
calculations for N = 3, . . . , 18 Pt layers (not shown in here) no significant N dependence of 
optical constants was found: the optical constants obtained for N = 3 are almost as good as 



those obtained with N = 18 Pt layers in the system. 

An extensive analysis of the inter- and intra-layer contributions to the permittivity e pq (u>) 
showed that in the case of periodic layered systems (simple lattices) the pq-like contributions 
to the permittivity are only functions of the relative position of layers p and q, i.e., e pq (u>) = 
i?\p-q\ (k,); the layer-resolved permittivities e p (to) are dominated by the corresponding intra- 
layer contribution e pp (to). Thus is not surprising that in the case of periodic layered systems, 
the following approximative form applies 



see Fig. 4. It was also found that the difference between the exact permittivity and the approx- 
imated one defined above, not only is independent of N, but also of uj. Thus a semi-infinite 
periodic layered system can be simply modeled as a sequence of seven layers (a central one plus 
three layers above and below this) in a constant dispersionless permittivity background, the 
contribution of which to the total permittivity is negligible. Because of the fast convergence in 
\p — q\ , calculations for periodic layered systems restricted to N = 3 layers can already yield 
reasonably accurate optical constants. 



It was shown analytically and numerically that in polar geometry and normal incidence 
the complex reflectivity coefficient of the right- and left-handed polarized light for periodic 
(layered ) systems is given by the Fresnel formula: only in the case of periodic (layered) systems 
the 2x2 matrix technique and the two-media approach become equivalent, i.e., they lead to 
identical Kerr spectra. This in turn clearly marks the limitations of the two-media approach, 
since applied to inhomogeneous layered systems, the corresponding Kerr spectra are dominated 
entirely by the optical activity of the surface layer. 

Inter- and intra-layer contributions to the optical conductivity for fee Ptjv/Pt layered sys- 
tems as calculated within the framework of the spin-polarized relativistic screened Korringa- 
Kohn-Rostoker (SKKR) method by means of a contour integration technique lead to optical 
constants, which are in very good agreement with available experimental data for photon ener- 




(18) 



Ip-?|=o 



SUMMARY 
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gies uj > 5 (5 being finite lifetime broadening/ The calculated optical constants for fee Ptjv/Pt 
do depend on the surface orientation; for N > 3 (number of Pt layers) they become virtually 
independent of N. 

The conclusion to be drawn for magnetic multilyer systems is rather simple: only the 2x2 
matrix technique that includes all multiple reflections and optical interferences can produce 
reliable spectra for this kind of systems. This in turn implies that a computational scheme 
such as the SKKR method designed especially for systems with only two-dimensional transla- 
tional symmetry has to be applied exclusively in order to evaluate the elements of the optical 
conductivity tensor. 
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FIGURES 




co (eV) 

FIG. 1. Reflectance of fee Pt bulk as calculated for the (100), (110) and (111) surface orientation 
(diamonds, circles and squares). The experimental data are taken from Refs. [20] (crosses), [21] 
(pluses), [22] (stars), [23] (triangle right), [24] (triangle left) and [25] (triangle up). 
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FIG. 2. Optical constants of fee Pt bulk as calculated for the (100), (110) and (111) surface 
orientation (diamonds, circles and squares). The experimental data of Weaver et al. [10] are displayed 
as grey circles. The dotted line marks the photon energy that equals the used lifetime broadening of 
5 = 0.653 eV. 
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FIG. 3. Permittivity of fee Pt bulk as calculated for the (100), (110) and (111) surface orientation 
(diamonds, circles and squares). The experimental data are taken from Refs. [20] (crosses), [21] 
(pluses), [22] (stars), [23] (triangle right) and [26] (triangle down). The dotted line marks the photon 
energy that equals the used lifetime broadening of 5 = 0.653 eV. 
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FIG. 4. Relative error made by approximating the permittivity of fee Ptjv|Pt(lll) by the sum 
over e r (uj) permittivities with r =| p — q |= 0, . . . ,3, see text. Up and down triangles, diamonds, 
circles, squares and stars denote data obtained for N = 3, 6, 9, 12, 15, 18. 
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